Nonlinear Current of Strongly Irradiated Quantum Hall Gas 
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Two dimensional electrons in weakly disordered high Landau levels are considered. The current- 
field response in the presence of a strong microwave field, is computed. The disordered Floquet 
evolution operator allows us to treat the short range disorder perturbatively, at any strength of 
electric fields. A simplifying Random Matrix Approximation reproduces the broadened Landau 
levels density of states and structure factor. We derive the magnitude of the Microwave Induced 
Resistivity Oscillations. The disorder short wavelength cut-off determines the non-linear electric 
fields of the Zero Resistance State and the Hall Induced Resistivity Oscillations. We discuss wider 
implications of our results on experiments and other theories. 

PACS numbers: 73.40.-c,73.50.Pz, 73.43.Qt,73.50.Fq 
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I. INTRODUCTION 

Some remarkable phenomena have been recently ob- 
served in high mobility GaAs/AlGaAs hetero-structures: 
Microwave Induced Resistance Oscillations (MIRO), 
Zero Resistance States (ZRS) 0, [H, 0, EL H, and Hall 
Induced Resistivity Oscillations (HIRO)[l, 0, @, 0, E3]. 
These experiments are carried out at weak fields (rel- 
ative to the Quantum Hall regime) and at tempera- 
tures T > hio c , (u> c is the cyclotron frequency), where 
Shubnikov-dc Haas oscillations are thermally smeared. 
Nevertheless, microwave radiation, or large Hall currents, 
can expose the underlying Landau quantization via the 
MIRO and HIRO oscillations. 

MIRO (in the Hall bar geometry) exhibits large mag- 
neto resistance oscillations as a function of radiation fre- 
quency u>, with nodes at harmonics of mu> c , m = 1,2,.... 
In strong enough radiation, the resistance in the positive 
detuning regimes may be nearly completely suppressed. 
This phenomenon is commonly denoted ZRS. The ZRS 
has been attributed to s pon taneous generation of inter- 
nal electric field domains [ll| . The macroscopic structure 
and stability of ZRS domains have been investigated by 
a phenomenological Lyapunov functional [l2l f]~3j ] . 

Microscopic theories for MIRO and ZRS, divide into 
two categories: (i) The Displacement Photocurrent (DP) 
mechanism, proposed by Ryzhii et. aZ.(l4| and Durst 
et. d.fjjfy. (See also Qj, G3 GJ|). Th e DP is stronger 
for well resolved Landau levels, whose width T is smaller 
than Landau level separation Huj c . (ii) The Distribution 
Function (DF) mechanism, proposed by Dmitriev et. al., 
[lj|. The DF mechanism becomes important when the 
inelastic lifetime is much longer than the elastic transport 
scattering time. A recent theory of HIRO for overlapping 
Landau levels was given in Ref. [20|. 

Independent of the particular dominant mechanism, 
the disorder potential is essential for the dissipative cur- 
rents. Treating disorder in high Landau levels is an old 
theoretical challenge: Since the clean system has macro- 
scopic Landau level degeneracies, a straightforward per- 
turbation theory in disorder is ill posed. The linear, 
dark conductivity, (in the absence of radiation), has been 



calculated by the Self Consistent Born Approximation 
(SCBA)[2l|,[23|. The SCBA is a selective diagram resum- 
mation, whose neglected vertex corrections are controlled 
by the smoothness of the disorder [23j|. Unfortunately, a 
controlled extension of the SCBA to strong static and 
time dependent fields, has so far proven difficult. Never- 
theless, the ZRS and the HIRO are inherently non linear 
effects which require theoretical attention. 

In this paper we devise a new 'divide and conquer' 
method [24J which incorporates strong electric fields with 
the disorder potential. It has long been appreciated that 
strong radiation effects on the clean Landau levels, are 
tractable using a Floquet transformation fLU [ItJ . How- 
ever, writing down an explicit Floquet transformation for 
a disordered Hamiltonian was believed to be intractable. 
Here we make progress by exploiting the commutation 
between Landau and guiding center operators. We con- 
struct a separable disordered Hamiltonian, which is com- 
pletely 'Floquet transformable'. This trick allows us to 
eliminate the electric fields, and obtain the random Flo- 
quet (quasi-stationary) eigenfunctions. These states, de- 
scribe our zeroth order densiity matrix, without any dis- 
sipative current. The dissipative current is subsequently 
computed perturbatively, to leading order in the remain- 
der short wavelength disorder. The physical small pa- 
rameter is the ratio of transport scattering rate to Lan- 
dau level width. 

This approach allows us to answer some outstanding 
questions: 

• What is optimal disorder for large MIRO and ZRS 
effects? 

We identify the 'figure of merit' as the ratio R = 
hu c /T, which is also expected to be large in systems 
exhibiting the HIRO effect. 

• What characteristic Hall fields determine the HIRO 
effect? 



The low intra-band field is 

E r = T/(el 2 B q s 



(1) 
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and the higher 'inter-band' field periodicity is given 
by 

E Wc = hw c /{el 2 B q s ). (2) 

q s is the high wave-vector cut-off of the disorder 
fluctuations, and Ib is the Landau length. 

• What is the expected magnitude of the spontaneous 
ZRS fields? 

The quantity Ep sets the overall scale of the ZRS 
field E zrs , which also depends on R, frequency and 
microwave power. E zrs has a second order dy- 
namical phase transition at a threshold microwave 
power. 

This paper is organized as follows: Section |TT] briefly re- 
views macroscopic transport theory of the Corbino and 
Hall bar geometries. The main purpose is to define the 
dissipative current, which is the subject of subsequent 
sections. Section IIIII introduces the microscopic model: 
Non interacting electrons in a weak magnetic feld, with 
'broad- tail' disorder correlations. The disorder poten- 
tial is split into two operators, where the long wave- 
length components are (mostly) incorporated into the 
broadened Landau levels Hamiltonian. In section IIV1 
the disordered Floquet evolution operator Ud{t) is con- 
structed. We show that while the Landau levels broaden 
into random Floquet eigenstates, the DC field (surpris- 
ingly) produces no dissipative current and a perfect clas- 
sical Hall current. The Random Matrix Approximation 
(RMA) captures the spectrum and eigenfunctions of the 
broadened Landau levels, and recovers the well known 
Self Consistent Born Approximation (SCBA). In section 
IVl we derive the leading order dissipative non linear cur- 
rent, Eq. P5|) . which can be disorder averaged numeri- 
cally. In section IVII the current formula is simplified by 
the RMA to a tractable analytical expression Eq. 
The simplified expression is analyzed in some detail: the 
SCBA dark conductivity is recovered, the MIRO, ZRS 
and HIRO effects are obtained. Predictions are obtained 
for the magnitudes of these effects and the values of non 
linear field scales. A plot of the full microwave irradi- 
ated current as a function of DC field is given in Fig. [5] 
We conclude with a brief discussion of theoretical and 
experimental issues pertinent to our results. 

Appendix [A] calculates the disorder matrix elements, 
and Appendix IBl derives the explicit Floquet operator. 



II. MACROSCOPIC TRANSPORT THEORY 

Here we briefly review magneto-transport theory of the 
Corbino and Hall bar geometries. This will provide a di- 
rect relation between the experimentally measurable cur- 
rents and voltages, and the quantity calculated theoreti- 
cally in later sections: the non linear dissipative current 
J d (E). 



In the presence of a microwave field E w , and a DC field 
E = BE, macroscopic transport theory assumes a local 
relation between the DC current density and the electric 
field, 

j = j d (E,E u ) E + j H (E), 
i H = (TflZxE. (3) 

j H is the non dissipative Hall current, and an is the Hall 
conductivity, an and j d (E) require input from a micro- 
scopic theory which involve quantum mechanical scatter- 
ing processes on length scales shorter than the inelastic 
dephasing length. 

In the presence of large scale inhomogeneities, (consid- 
ered in Refs. [H, G3 ) j the dissipative current j d (r) is not 
necessarily parallel to E(r), and an may vary in space. 
Here we do not consider such large scale disorder. 

A. Corbino Geometry 

By the definition ([3]), for a position independent an, 
the dissipative and Hall currents are conserved sepa- 
rately: 

V-j = V-j H = V- (/E) =0. (4) 

In a Corbino geometry (see FigO}, one specifies the po- 
tential on the inner and outer boundaries of the sam- 
ple which determines the purely radial electric field E = 
E(r)r. The dissipative current j rf is also radial, while the 
non dissipative Hall current circulates around the annu- 
lus unrestricted. 

The Corbino geometry completely separates the effects 
of (jff and j d , and is simpler to analyze theoretically. 
The function j d (E) is uniquely determined for particular 
sample dimensions, from the total current and voltage 
between the outer and inner edges, with no dependence 
on o~h- 

As was shown in Refs. fill. Il2l [l3j . in the absence of 
an external bias current, a negative conductivity in some 
field range j (E) < may be unstable to a spontaneously 
created domains of finite internal fields whose magni- 
tudes are |E| = E zrs and orientation is determined by 
the boundary conditions and long range potentials. 

f{E zrs ) = 0. (5) 

In the absence of 'pinning' in the Lyapunov functional, 
the domian walls are free to move and absorb any charge 
in external bias voltage leaving the current vanishingly 
small. This is the so-called 'Zero Resistance State' 
(ZRS), which, confusingly, yields zero conductance in the 
Corbino geometry. 

B. Hall Bar Geometry 

Experimentally, the Hall bar geometry (see FigfTJ) is 
more popular for ease of fabrication. The Hall bar ge- 
ometry imposes j y = 0. The current is forced through 
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FIG. 1: Corbino and Hall Bar Geometries. In the homoge- 
neous Corbino sample, the electric field and dissipative cur- 
rent density are radial, while the non dissipative Hall curent 
is azimuthal. In the Hall bar, the current between leads in- 
cludes both dissipative and non dissipative components. For 
large Hall angle oh » cr X x, most (but not all) of the current 
through the bar is non dissipative. 



in one direction j = jx, and the longitudinal and trans- 
verse voltage drops are measured. The electrochemical 
field E = (E x ,E y ) can be deduced from the geometry. 
The can be deduced. By ((3]), the electric field satisfies 
the coupled non linear equations 

E x 



3 = 3 d (E)^+a H E y , 
= -a H E x +j d (E)^-. 

In the case of large Hall angle, 

o H » o xx = j d /E, 



(6) 



(7) 



and thus E m E y » E x . This inequality simplifies the 
solution of ([6]): 



Ey(j) 



J d (E) 

OH 



(l + 0{\a xx \/a H )) 



(8) 



where we can drop the terms of order (j xx /<jh << 1- 
An experimental measurement of E x (j) approximately 
describes the non linear function j d (E), which will be 
calculated in Eq. |48)) . 



(9) 



eaj) « —r (-*- 

oh \oh 



The negative conductivity implies E • j < 0. In the ZRS, 
the spontaneous fields E zrs of ([5]), arc associated with 
disssipationless flow of Hall currents with magnitude 



oh E z 



(10) 



In the ZRS, the currents traverse the system in the x 
direction producing very little longitudinal voltage drop 
and dissipation. For the Hall bar geometry of width W 
(see Fig. [T]), the resistivity remains vanishingly small 
upto a critical current I cr given by 



Wa H E z 



(11) 



Thus, the critical current directly measures the internal 
ZRS field. 



III. THE MICROSCOPIC MODEL 

We consider a two dimensional electron gas (2DEG) 
subject to a perpendicular magnetic field Bi, microwave 
field E w , and a DC field E^ c . The standard treatment of 
electrons in a magnetic field defines dimensionlcss Lan- 
dau level and guiding center coordinates as 



Ib 



1 



-z x r, 



R 



r 

21b~ 



2l B 
— z x p, 



(12) 



where the magnetic length is Ib — y/hc/eB. By con- 
struction, the two guiding center components commute 
with the two Landau operators: 



[n a ,R ]=O, a,P = x,y. 
The microscopic model Hamiltonian is 



H(t) = 



a(*)r + V(r), 



(13) 



(14) 



where uj c = eB/(mc) is the cyclotron frequency. The 
electric fields are represented by the dimensionlcss gauge 
field 



a(t) 



h 



(ft 


jjjg — hat 











(15) 



The disorder potential is described by a random function 
of position r, 



q 



(16) 



where A is the area of the system. The Fourier compo- 
nents are random complex numbers obeying Vq = V*^, 
with correlations 



(VqV- q ,)=AW{q)8, hq > 



(17) 



Our 'divide and conquer' tactic is to split the disorder 
potential into two components by a dividing wavevector 
qi. We assume that W(q) is a 'broad-tail' distribution, 
(see Fig. [5]), i.e., most of its weight is concentrated in 
a region below some dividing wavevector qi, (say, of the 
order of Ib), but the higher moments are dominated by 
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much shorter wavelengths. For example, we might choose 
the function 



W(q) = Wie- 2 */* +W a 0(q a - q), 




m q Qs 

FIG. 2: The 'Broad- Tail' model of disorder correlations. The 
wavevector qi bounds the long wavelength spectrum which is 
incorporated into the disordered Floquet evolution operator 
(see text). Most of the spectral weight is below qi, although 
the higher moments which determine the dissipative current, 
are dominated by the short wavelength cut-off q a - 



This model is realistic for high mobility heterostruc- 
tures, where the offset distance to unscreened charges is 
d = 1/qi [25|, and W s is the intrinsic short range disor- 
der of the 2DEG. The upper cut-off is effectively at some 
q s < ixkp. We do not include at this level any electron 
phonon or electron electron interactions. 



where m = 1, . . . til- The state \k) is centered at position 
(y) = klg. The matrix elements of V; are the (identical) 
block-diagonal matrices (see Appendix |A|): 

V& = \ E V^L nF (£) e"-^<W, 

A |q|<* V J 

(21) 

L n is the Laguerre polynomial of order n, and q = qlp- In 
pop . e a and tp a (k) = (k\na) are the disordered spectrum 
and eigenvectors respectively. 

V;, by construction, has no inter-band matrix elements, 
and identical intra-band blocks. It is therefore a random 
function of only the guiding center operators R (see Eq. 
(fT2"]) . Nevertheless, it contains most of the spectral weight 
of the disorder. 




kl B 

FIG. 3: A typical random eigenstate ip a , chosen near the 
broadened Landau band center. The phase (amplitude) is 
plotted in the top (bottom) panel. Mb describes the height 
position on the cyclinder. We use Landau degeneracy np = 
128, Fermi level rip = 15, and cut-off qi = 2/Ib for the long 
range disorder. The cyclotron radius and Landau length are 
depicted. 



A. Splitting the Disorder 



The main trick is to split the disorder into two opera- 
tors 

V(r) = V,(R)+V s (tt,R) (19) 

The long wavelength disorder term is defined as the op- 
erator 

Vi = V^\nk)(nk'\ =J2^\na)(na\. (20) 

nkk' na 

up is the single-spin filling factor at the Fermi energy. 
The Landau level degeneracy is np = A/ (27rZg), where A 
is the system area. Without loss of generality we choose 
the basis \nk) on a cylinder of circumference L x = y/A, 
where n is the Landau level index, and k — 2itmlB/ L x , 



B. The Random Matrix Approximation 

We shall now see that at high Landau levels, V£j£ 
is well represented by a random Hermitian matrix. Its 
eigenstates appear as random extended wavefunctions as 
depicted in Fig. [3] 

Although numerical averaging over disorder poses no 
great difficulty, the statistical properties of V£j5 suggest 
to approximate it as a Hermitian random matrix of the 
Gaussian Unitary Ensemble (GUE). This Random Ma- 
trix Approximation (RMA) will provide easy averaging 
of spectra and correlations over the long wavelength dis- 
order. This presents another example of the usefulness of 
Random Matrix Theory [27[ for quantum transport [2£| . 
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Fortunately, the RMA becomes exact in the limit of 
a uniform white-noise disorder fluctuations W(q) = W, 
and large filling factor. The matrix elements of V k "£ be- 
come maximally uncorrelated: 



ttR c A 

x 

/2ttZ 2 



—y 



W(q) 

I 



Mx(k-k"') 



q 



X <5fc',fe+q H ^fe",fe'"+g 



f 



RMA 



2ir 2 R, 



dqW(q), 



(22) 




r/r 



.RAM 




2T* ,1 



(rtj? — l)7iijj c 



Qih 




c JO 



where R c = lB\J2np. The density of states is defined as 
an ensemble average 



>{E) = (£5{E -E na )) 



Vr 



(23) 



FIG. 4: Numerical averaging of 500 disorder potentials, for 
the density of states p(e) of Vuy , using Wi — 1%, uf = 80, 
Landau level degeneracy til — 512 and qi = 2/Ib- The fit is 
to a Gaussian of deviation F Inset: dots are the numerical 
values of F for different cut-off wavevectors qi , divided by the 
Random Matrix Approximation r RMA , Eq. ([22)) . 



For a GUE ensemble, the density of states of each Lan- 
dau band is expected to be semi-elliptical. In the limit 
of narrow Landau levels, the SCBA also produces semi- 
elliptical density of states. Moreover, the SCBA [HI, [23| 
for the Landau level broadening (sometimes called in 
the literature 'quasiparticle scattering rate' h/r s ) equals 
Frma, within the assumptions of Eq. (fT5|) . 

Here, however, we are interested in long range disorder 
where qi << kp < oo. This implies that the matrix 
elements of V£ff are correlated within a distance Ak < 
1/(IbQi), and that the density of states at the tails will 
not be well described by the RMA's semi-elliptical form. 

In the regime qils € [0.3,5] the numerical density of 
states, averaged over a Gaussian distribution of disorder 
potentials {F q }, can be fitted to a sum of Gaussians: 



HE) 



1 



1—^ (E-nhu, a )* 

e ar 2 



(24) 



As shown in Fig|4]the Gaussian fit is excellent for qils = 
2. The fitted Landau level width T is also quite close to 
the RMA result, 



r « i.3r 



It MA 



(25) 



The RMA for the eigenfunction correlations, which de- 
termine the structure factor, will be investigated later in 
Section [VDj 



IV. FLOQUET TRANSFORMATION OF TIME 
DEPENDENT FIELDS 



A. Clean Landau Levels 



The clean Hamiltonian, driven by a time dependent 
vector potential is 



W)~k + a(i)| 2 



(26) 



The Schrodingcr equation for the evolution operator 
ihd t U F (t,to)=Ho(t)U F (t,t ), (27) 
is solved in Appendix [B] and is given by 



- p -^(t)-7T p -ifio(t-to) p +i£(t )-TI 



U F (t,t ) = e 

Ho(t) = hnuj c \nk) (nk\ 



(28) 



where the Floquet Hamiltonian 7Yo is identical to the 
time independent Hamiltonian. This identity, (called 
Kohn's theorem) also holds in the presence of arbitrary 
two-body interactions J2i<j ^( r « — r j) between electrons 

The vector fields £(i) = £^+£0 arc solutions of coupled 
linear equations. Their explicit form, given in Appendix 
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iBl is well known [il 



ftio 



3t 



w c E», + iuE, 



Huj 



Ut) = t^V dc xz-^-E dc . (29) 
h to. 



B. Floquet Transformation With Disorder 

In Eq. (|2"D|) we have constructed the long range disor- 
der operator V/ to a function only of R. Thus, Eq. (flU]) 
implies the trivial, yet very useful, commutation relation 



[Vi,n}=0, 



(30) 



which we exploit to write the time dependent disordered 
Hamiltonian as 



n d {t)=Ho(n,a(t)) + Vi[TL} 



(31) 



Now we can construct the disordered Floquet evolution 
operator U d , which obeys the evolution equation 



ihd t U d (t,t ) =Hd{t)U d (t,to). 



(32) 



Since the Floquet operators e z ^')' 7r commute with V; [R] , 
the solution is simply 

U d (t,t ) = e-^W-w e -iW*C*-*o) e +<€(to)-W j 

i-i d — (hnuj c + e a ) \na) (na\ , (33) 

na 

with the same vector function £(£) as for the clean Lan- 
dau level problem (|2"9"|) ! For Ed c =0, the Floquet states 
\na(t)) = U d (t,to)\n,a) return to themselves (up to a 
phase) at every integer period. 

Expression (f3"3")) implies that the Floquet quasi-encrgics 
(the spectrum of the transformed Hamiltonian TL d ) is the 
disorder broadened Landau levels. The equivalence be- 
tween stationary and quasi-stationary states with time 
dependent electric fields, is a special property shared by 
harmonic oscillators, Landau levels, and our H. d . 

Using the last equation in (|B3|) , we can also explicitly 
express the time dependent Fourier operator as 



e -K n d t e i ^ T e - J k' H di ft <Jsq-£(t) 



(34) 



where we can use the explicit solution (|29[) to perform 
the Bessel expansion, 

oo 

e «sq-£(t) _ e -iV jy j v (qA U q) e if ( I/w +( ei B/ R )q xE rf<:- 2 ) ; 

(35) 

where "0 is an unimportant phase, and J v is the Bessel 
function of order v. The quantity 



A wq (E w ) = max t [q •£(£)] 



(36) 



is proportional to the microwave field strength. For pos- 
itive circular polarized light of field E w = (x + iy) , £ 
and A^q diverge at the cyclotron resonance as \u— w c | _1 - 
For the opposite polarization, A^q is continuous at res- 
onance. 



V. CALCULATION OF THE CURRENT 

The DC current density is given by the full evolution 
operator U as 



eh 



Amli 



J2pna(na\W(t,t )iTU(t,to)+a(t)\i 



(37) 

where p na is the electron density operator, defined in the 
Floquet basis \n, a). 

The notation F(t, to) implies averaging the end-point 
times t, to over many periods of 2tt/uj. 



A. Zeroth Order Theory 

The zeroth order DC current density is given by setting 
U -> Ud in |37|) . Using (|B2j) and l[29|). it is easy to obtain 



j<°> = (| xz + a dc ) 




mis A 



e 2 n e 
muj c 



z x E 



d c ? 



(38) 



where the result holds for any distribution p na (not nec- 
essarily Fermi-Dirac) which sums up to the full electron 
density n e . 

It is interesting that although H d has random spectrum 
and wavefunctions, its current density (|38p produces no 
dissipation. In addition we see that the Hall conductivity, 
is precisely given by the classical value, 



0~H = 



2 2 

e n e e 
= nF-r 

mu) c h 



(39) 



although for V/ ^ 0, our system is by no means Galilean 
invariant. Incidentally, this result (|38p applies in the 
presence of arbitrary electron-electron interactions which 
may be added to Ttd(t). 

We iterate that Eq. (f38|) applies only to the zeroth 
order Hamiltonian defined in Eq. (f3"Tj) . Dissipative effects 
enter at second order 0(Vg). 



B. Second Order Dissipative Current 

The leading order dc current, requires expanding the 
momentum operator to second order in V s . 
The interaction representation is defined as, 



O d (t, t Q ) = ul(t,t )OU d (t, t ). 



(40) 
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The zeroth order ('clean limit') Landau operator (in 
complex notation 7r = tv x + i-n v ) transforms as 



(41) 



Averaging (|43|) over {U q } keeping V; fixed, yields 

(„oi^.( t , t0 ,w,^i:*±* 

qm 



where 7 is an unimportant c-number function. 

The full evolution operator U is expanded to second 
order in V s , 

U(t,t ) = U d (t,to) I *' V «(*''*o) 

,H) 2 '* 



W(q) (1 - S ) Inm(t, to), 



Inm(t,t )=J2 / dt ' 



t-to 



Jui c (t-t') 



±i(E na -E mot ,)(t' -t") 



K 2 



to J to 



dt' / dt"V d (t,t )Vf(t",t ) + 



(42) 



(47) 



which yields the second order correction to the Landau 
operators as 

t t' 

Tt^(t,t ) = -=, I dt' [ dt"e iu ^-^ 

h Jo Jo 

x [[*M d) (t',to)] M d \t"M 

t t' 



r (d) (t',*o) e -»q'-r (d >(t",to) 



(43) 



The leading order dissipative current is obtained by 
the disorder averaged expression, 



Ami 



B 



^2fna ((?ia\TvW(t,t )\na] 



Here we use the Fermi-Dirac distribution 

-1 



(44) 



(45) 



Corrections to the distribution function enter at higher 
order in V s , see discussion in Section fVIIl item 2. Eq. 
(|44|) ignores correlations between Vi and V s , which exist 
in the range q < qi. The resulting error in the current, 
which is dominated by large wavevectors, is relatively 
suppressed by powers of qf /q 2 . 

Using the matrix elements of the time dependent 
Fourier operator (f3~4"|) we obtain, 

{na\e^- rW ^\m0) = e i( E --*W)(*-*°) £> nm (q) 



k 



^2 Jv (^ A ^q) 1 



y— — 00 



(46) 



We can henceforth neglect the terms (—S n mO q < qi ) in (|47|) . 
The corrections contribute only to the low wave vector 
intcrgration, and are therefore suppressed by factors of 
(«/?«) « 1- 

Integrating (j4"4")l over time, we obtain the the full cur- 
rent expression, 



mR c io c 



d 2 q 



(z x q) 



W res (q) 



J dede' (/„, e - f m>e ,) K 2 (q, e, e') £ | J„(gA wq ) 



x5{E me > - £; ne - 7u>w - e^|z x q ■ E dc ). 



(48) 



is the disorder averaged intra-band structure factor, 



K 2 (q,e,e') 

1 y 



dtdt' 
{2nh) 



Jqx(k-k') 



,-e» 



.4 



x (<Pl(k)<Pa(k')<pp(k + q y )(pp{k' + q y ) 
x 5(e-e a )6(e' -ep)) v 



(49) 



Note that by the symmetry of / d 2 q integration, the sec- 
ond order current is parallel to the DC field, and purely 
dissipative. Hence to second order, the Hall conductivity 
remains equal to the classical value o~h of (|39[) . 



C. Simplification at T > hu> c 

Many of the relevant experiments are carried out at 
temperatures T > ftw c , where the Shubnikov-dc Hass 
oscillations are thermally smeared. Using the energy- 
conserving delta function, we approximate the difference 
of Fermi functions as 

/n,e - f m>e > « -0 (T - \AE V \) AE v df ^ € \ (50) 
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where AE U is the transition energy 

/S.E V = hvuj + el%z x q • E dc . (51) 
Dchning the single frequency structure factor K\ as 



Ki (g, u>) = J deK 2 {q,e,e + hw), 



(52) 



the current expression simplifies to 



/ d 2 q — ^(z x q) 

J q 



\AE V \<T) 

x {ftvu) + el%z x q-E dc ) |J„(qA W q)| 5 



(53) 




one has 



-l(^5( e -e Q )5(e-e a )) = 27rZ|p(e)p(e') 



a/3 



(54) 



for e — e' >> The RMA approximates the dynamical 
structure factor Eq. (|4"5)) as 



A". 



KM A 



1 cxpf-^), (55) 



(27r)»i|ra 



2r 2 



where T = r flMA is the width of the single Landau level 
density of states (|2"2"|) . 

The RMA for the single frequency structure factor is 



deK? MA 



(e, e + 



exp 



4r 2 



(56) 



We have computed the structure factor K\(q, u>) numer- 
ically for a finite Landau level degeneracy on a cylinder. 
As depicted in Figs. [5]and[6l the RMA works well for 
q >> qi, where T is given by (|25[) . This finding con- 
firms our expectation that the RMA is applicable at large 
wavevectors, since correlations between matrix elements 
of V^J, decay beyond the wavevector scale qi. In con- 
trast, at low wavevectors q < 3qi, there is a peak in K\, 
which is absent in the RMA. It contributes negligibly 
however, to the overall current. We note that there is 
no signature of enhanced scattering at wavevectors cor- 
responding to the classsical cyclotron diamater 2R C , see 
discussion in section fVIIi item 6. 



FIG. 5: The static structure factor K\{q,u> = 0) evaluated 
numerically for ul = 128, uf = 15, and qi = 1/1%. The 
exact correlator has a peak at zero waavevector of width qi, 
and agrees with the (wavevector independent) RMA at higher 
wavevectors. The Inset depicts classical cyclotron orbits, 
which strongly overlap when separated by 2R C . In contrast, 
K\ shows no signature of enhanced scattering at wavevectors 
corresponding to 2R C . 



D. The RMA For The Structure Factors 

Within the RMA, the Hamiltonian of each Landau 
level is taken as a member of the Gaussian Unitary En- 
semble (GUE). Thus, by the RMA, the wavefunctions 
(p a (k) are normalized complex random vectors, and the 
spectral correlations are completely uncorrelated beyond 
a few level spacings |27| . From random matrix theory, 



VI. ANALYTICAL RESULTS 

The full non linear current j (E^o E w , u>, B, T, n e ) can 
be computed from Eq. (|48p . by numerically averaging 
K\ over disorder. It is useful however to use the RMA 
and obtain analytical results. These allow a better un- 
derstanding of the important characteristic energy and 
field scales of the non linear current. 

The first order of business, is to extract the dark linear 
conductivity from Eqs. ([53"]) and (j56"|) . We obtain: 



Hark 



ripe 



2ir2T tr T 



(57) 



Not surprisingly, this result agrees with the SCBA [291 ] . 
where the transport scattering rate is defined in the stan- 
dard way : 



1/rtr = {21th 2 VFkp) 



dqq 2 W(q). 



(58) 
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FIG. 6: The dynamical structure factor. Kf am {q, uj) is eval- 
uated numerically, for til = 64, uf = 80 and qi — 1/Ib- The 
Random Matrix Approximation K^ MA is given in Eq. (|56[l . 
Good agreement is found for q > 3qi, at all frequencies. 



Using the RMA for K\ , given by (|56|) , the current can 
be simplified into an analytical expression: 

j d = a^ lk E r (^ dc F(e dc ,e u ,u;) 

+ (^jG(s dc ,e u ,uyj. (59) 

The dimensionlcss fields are defined as, 
£dc = E/Er, 

£ w = A W q,9 s / B . (60) 

The characteristic electric field scale 

E r = T/(el%q s ), (61) 

defines the lowest scale of non-linearity in j d . 

At |E| « Er, resonant intraband transitions are en- 
abled over a distance of l%q 8 - 

For a white noise short range disorder W(q) = 
W s & q< q s , the dimensionlcss functions F, G in (|59|) are ex- 
plicitly evaluated as 

|AB„|<T 

F = ~ E / d^Ure^ 2 

x e -(ye dc +h(vuj-m.LU c )/r) 2 /4 

|A£„|<T 

G = — V 1/ / d 2 r^|J„(r £w )| 2 

77 — ' Jr<l r 

7TLV 

x g-(l/e£ic + 'i(i / w-™"c)/r) 2 /4 

(62) 



where the sums arc limited by the temperature condition 
on the transition energies (|51D . 




-4-2 2 4 

oj — ifioj c [r/?i] 



FIG. 7: MIRO effect. The frequency dependent conductivity 
is plotted against detuning frequency near the m th resonance. 
A set of dimensionless microwave field strengths are given. 
The broadening ratio (|65jl. is chosen at R = 10. Tempera- 
ture is T = 7hu> c . The regions of negative conductivity are 
unstable toward the Zero Resistance State. 



A. MIRO: Photo-conductivity Oscillations 

G describes the photo-current contributions. Near a 
particular primary [y = 1) resonance at some interband 
interval m > 1, the dimensionless detuning frequency is 
defined as 

<L = ^(w - fhu c ). (63) 

The MIRO effect, arising from the function G, changes 
sign at the resonances as 

G oc -<L- (64) 

The magnitude of MIRO is controlled by the 'figure of 
merit', the Landau level broadening ratio 

T V Wiqimc J 

For R >> 1 the total current at 5 U > 0, can turn neg- 
ative and produce the Zero Resistance State. This ratio 
increases as \f~B~. Our approximations, however, break 
down at low filling factors. 

In Fig. [7J we plot a low field conductivity oscillation 
as a function of detuning frequency. We choose R = 10 
as the broadening ratio, and a set of microwave radiation 
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field strengths e u . At temperature T = 7frui c , we note 
that the conductivity at 5 U = 0, is independent of ra- 
diation power. The regions of negative conductivity are 
unstable to formation of a ZRS. 



B. HIRO: The Dark Non-Linear Current 

In the absence of radiation, G vanishes, and the non 
linear conductivity in (|59|) is <r^ r k F. In Fig. HI we plot 
F{E) for different values of hui c /T. 

At weak fields of order E-p , we see a Gaussian decrease 
of conductivity due to diminishing phase space for intra- 
Landau bands scattering. Caution has to be exercised in 
trying to fit the low conductivity regime with Eq. (|48|) . 
The far tails of the density of states require multiple scat- 




10 20 30 

E [Br] 



FIG. 8: The HIRO effect. Non linear current oscillations at 
zero radiation, as a function of at DC electric field, for a set of 
broadening ratios R. The field scale is Er is defined in (|61|l . 
The oscillations periodicity is determined by E Uo = R Er- 




10 20 

E [E T ] 

FIG. 9: Non linear current-field characteristics for differrent 
levels of microwave field e w . The field scale is _Er l|6ip . Tem- 
perature is chosen at T = 7fiui c , and the broadening ratio is 
R = 10. Su,, defined in (|63|l . is fixed at 0.9. Negative current 
regimes are unstable to a ZRS with spontaneous electric fields 
of magnitude E zrs . 



4, =0.9 

_ m • 
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FIG. 10: The ZRS spontaneous field. E~ rs is given by the 
zero of j d (E) as a function of dimensionless microwave field 
e W) at fixed detuning frequency 5^. A continuous (second 
order) dynamical transition occurs at a threshold microwave 
field. 



At larger fields the conductivity has secondary maxima 

at 

|E maa; | = (m + S)E LUc , (Jps 0.25, m= 1,2,... 
E^ c = -^=RE r , (66) 

which corresponds to inter-Landau level scattering over 
the length scale q s l%- For q s > nkp, by (|A4[) . one must 
use q s — > kpir. 



C. The Non Linear Photocurrent and ZRS Fields 

In Fig. [9l we plot the fully irradiated non linear cur- 
rent at different microwave fields, at one fixed optimal 
detuning of 8 U = 0.9. For weak radiation, the HIRO os- 
cillations are observed, and their phase shifts as the mi- 
crowave power increases. At stronger radiation, the neg- 
ative conductivity regions are unstable toward the ZRS. 
The ZRS creates spontaneous electric fields of magnitude 
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E zrs , where the dissipative current vanishes j d (E zrs ) = 
[13, El. In Fig. [TUl wc plot the dependence of E* rs {eJ) 
at fixed detuning frequency. The ZRS exhibits a con- 
tinuous dynamical phase transition as a function of mi- 
corwave power, which was previously predicted on phe- 
nomenological grounds, by Alicea et al. [3l| . 



VII. DISCUSSION 

We emphasize several issues regarding our results. 

1 . Our analysis is limited to second order in the short 
range disorder V s , which is controlled by the small- 
ness of W s q s /(Wiqi) < 1. In physical terms, this 
is equivalent to the well known 'small angle scat- 
tering' limit where h/(Tt r T) < 1, with T and Tt r 
defined by Eqs. ([25]) and ([58]) respectively. Small 
angle scattering is used to justify the SCBA [23| . 
By Eqs. ([55)1 and (|57|) . the small angle scattering 
parameter can be experimentally determined by the 
Hall angle 6 H , 



cotan(0#) 



T dark 



Off 



27T2 T tr T 



(67) 



By Eq. ([7]), a large Hall angle 9h ~ tt/2 is also 
the limit which effectively separates the effects of 
dissipative and Hall fields in the Hall bar geometry 



dissipe 



2. Our photo-conductivity supports the Displacement 
Photocurrent mechanism [3, EH- We have as- 
sumed here, that electron-electron interactions are 
relatively weak for well resolved high Landau levels. 
However, we also implicitly assume that electron 
phonon interactions are strong enough to produce 
an inelastic scattering time of the order 



Ttr 



(68) 



Hence, to the leading order in l/r tr , we do not need 
to correct the Fermi Dirac distribution. We do not 
rule out that DF mechanism [n| may yield sizeable 
corrections, for the case Ti„ > r tr . 

3. We have assumed a model of weak disorder which 
produces well separated Landau levels. We find 
that a large ratio hu> c /T, can produce large MIRO, 
and ZRS effects under radiation. It also is respon- 
sible for a large HIRO effect. Indeed these effects 
have been observed in the same samples [Toj . 

4. The conductivity at lo = muj c appears to be insen- 
sitive to the radiation field (see Fig. [7]). Within Eq. 
(|59[) . this is true only at temperatures T/(hw c ) = 
v max >> otherwise the Bessel functions sum 
rule, 



lim 



E 

\u\<V™ 



\J„(re^)\ 2 = 1, 



(69) 



is not saturated. At very low temperatures, only 
a few Bessel functions contribute, and we expect 
the conductivity at resonance frequencies to be sup- 
pressed under radiation. 

5. The characteristic field E-p (|6"T1) . is in the same ball 
park as Ryzhii's estimate [Tj]- It differs from the 
much larger estimate (> Hu) c / (eR c )) provided by 
Ref. [H. 



6. The characteristic field of secondary HIRO max- 
ima ([51]) depends on the disorder cut-off q s . This 
is in contrast to the semiclasssical picture @, [2(| 
of 2kp scattering, or touching cyclotron orbits of 
radii R c = vf/lo c . This scenario is depicted in the 
inset of Fig. O This argument suggests enhanced 
scattering at an electric field 



E- 



2R C 



hxo c 
2eR, 



7^ Eu. 



(70) 



where our E Uc in (|66[) . depends on q s . As evidenced 
in Figs. [3l and [5] there is no signature of enhanced 
scattering at wavevectors 2R c /l^, unless, by coin- 
cidence, the short range cut-off precisely matches 
q s = 2k p. In summary, the cyclotron radius is not 
a noticeable lengthscale in the structure factor at 
high Landau levels. 

Experimentally [H, values of oscillation field have 
not been universal, and have varied in the range 
(1.63,2.18) x E 2Ra . 

7. An important issue, which we have ignored in this 
work, is the relation between the dimensionless mi- 
crowave field e u and an external microwave radia- 
tion field and polarization. The resonance of pos- 
itive circular polarization at lo = oj c implies that 
close to the cyclotron resonance, one cannot ignore 
strong frequency and polarization dependence of 
the dielectric function. There are related open is- 
sues raised by recent experiments, most notably the 
apparent independence of the MIRO on microwave 
polarization [32j, and the seeming inability of 2 pho- 
ton absorption to explain a MIRO about uj c /2 j33l |. 
We shall defer discussion of these problems to fur- 
ther investigations. 
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APPENDIX A: THE DISORDER MATRIX 
ELEMENTS 



Using the definition (|12|) . the commutation between 
Landau operators and guiding center coordinates are 
readily obtained, for a, ft = x,y 



[n a ,Rp\ = 0, 

[7Txi%] = [Ry,Rx] = i- 

The full disorder operator is expressed as 



(Al) 



(A2) 



by 



The matrix elements of the Fourier operators are given 



(n\e 



= e — a — IqZ^/V^I 



-•iisq-Vrxz I 



(m— n) i{m—n) 



xW-L™-"((^) 2 /2) 
V n! 

£>nm(q), 



(A3) 



where </> = arctan^/tfc). Eq. (|A3j) is used to split the 
long and short wavelength terms in V in Eqs. (|21[) and 



The asymptotic expression for the associated Laguerre 
polynomials [34j valid at large n < m, with q = qlg < 
v2jwt, yields 



9pnm(r/2)| 



2n7r 



cos 2 (V2ng - <t> n - m ), (A4) 



where </> = 7r/4 + 7r(m — n)/2. The fast oscillations of 
the cos 2 can be replaced by its constant average of i, 
which simplifies the q integration over smooth functions 
considerably 



APPENDIX B: THE EXPLICIT FLOQUET 
OPERATOR 



The full disordered evolution operator (|28|) is decom- 
posed into 

U F (t,t ) = W{t)exp(-^Ho{t-t )) W*(t ), (Bl) 

where TIq is the clean Landau level Hamiltonian. and 
the vector-potential-dependent Floquet operator W is 
defined as 



wit) = e -^( t )- 7r e - i / tL (*') d *'. 



(B2) 



Here we shall determine the fields £(t),L(t). W acts 
on the time derivative, and translates the Landau level 



coordinates 7r as follows 

W f (id t ) W = 

W^IVW = 7T + |XZ 

W f rW = r + Z B £, 
which transforms the evolution equation as 



(B3) 



id t 



77 + ^xz + a- £/u c 



\k\ 2 1 

= ihd t - Hd- 



(B4) 



The fields £ are given by solving the linear equations 



£ - u c £ x z = w c a(i), 



(B5) 



For convenience, we fix the origin of time t = 0, such that 

€(* = o) = €o- 

Using the complex notation £ = £ x + i£ y , we write 



£ + iu; c £ = LO c a(t). 



(B6) 



The dimensionless gauge field (written in circular polar- 
ized components) is, 



Qj'X, I XCly 

if: 

h 



el B (iE+ 



Separating the solution into = £ + + £ 
obtains 



• (B7) 
fo, one 



±- 



elnLO r E 



.elsEdc elsEdc 



-t 



ftlJr 



(B8) 



The zero frequency component of £ in cartesian coordi- 
nates is 

€o(*) = ^ Vdc x z - -^E dc . (B9) 

In vector notation, we thus derived the explicit form 
of the Floquet fields £(£) as given in the main text by 
Eq. ~ 
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